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ABSTRACT 

Recent measurements by the Pierre Auger Observatory suggest that the composition of 
ultra-high energy cosmic rays (UHECRs) becomes dominated by heavy nuclei at high ener- 
gies. However, until now there has been no astrophysical motivation for considering a source 
highly enriched in heavy elements. Here we demonstrate that the outflows from Gamma-Ray 
Bursts (GRBs) may indeed be composed primarily of nuclei with masses A ~ 40 - 200, which 
are synthesized as hot material expands away from the central engine. In particular, if the 
jet is magnetically-dominated (rather than a thermally-driven fireball) its low entropy enables 
heavy elements to form efficiently. Adopting the millisecond proto-magnetar model for the 
GRB central engine, we show that heavy nuclei are both synthesized in proto-magnetar winds 
and can in principle be accelerated to energies > 10 20 eV in the shocks or regions of mag- 
netic reconnection that are responsible for powering the GRB. Similar results may apply to 
accretion-powered GRB models if the jet originates from a magnetized disk wind. Depending 
on the precise distribution of nuclei synthesized, we predict that the average primary mass 
may continue to increase beyond Fe group elements at the highest energies, possibly reaching 
the A « 90 (Zirconium), A « 130 (Tellurium), or even A « 195 (Platinum) peaks. Future 
measurements of the UHECR composition at energies > 10 20 eV can thus confirm or con- 
strain our model and, potentially, probe the nature of GRB outflows. The longer attenuation 
length of ultra-heavy nuclei through the extragalactic background light greatly expands the 
volume of accesible sources and alleviates the energetic constraints on GRBs as the source of 
UHECRs. 
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1 INTRODUCTION 

The origin of Ultra-High Energy Cosmic Rays (UHECRs) 
is o ne of the great mysteries in high energy astro physics 
(e.g. iBlandford & EichlejTl987l : iNagano & Watsoidl2000h . UHE- 
CRs are generally thought to originate from extra-galactic dis- 
tances: the break observed in the cosmic ray spectrum at ~ 3 x 10 18 
eV (the 'ankle') is often interpreted as the energy beyond which 
the Galactic magnetic field can neither isotropize UHECRs nor ap- 
precia bly prolong their residence time in the Galaxy (e.g. iHillasI 
2005). An extra-galactic origin is also suggested by the cut-off 
observed in the spectrum above ~ 6 x 10" eV (Abraha m et al.l 
2008; lAbbasietal .1 12008). which is generally interpreted as the 
result of UHECRs (protons or heavy nuclei) interacting with the 
cosmic microwave background (CMB) and other sources of extra- 
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galactic background ligh t (EBL). This is the ' GZK' effect, initially 
propo sed for protons by iGreisenl l i 19661) and IZatsepin & Kuz'rninl 
dl966h . 



Only a handful of astrophysical sources are plausible sites for 
accelerating UHECRs because the requirements on the magnetic 
field, compactness, and energy budget are stringent. Commonly 
discussed candidates can be divided into persistent and transient 
sources. Persistent sources include powerful relativistic jets from 
Active Galacti Nuclei ( AGN: e.g. iMannheim & Biermannlll992l: 



Berezinskv et al.l 120021: iFarrar & Gruzinovl |2009| ; iDermer et alj 
20091: lTakami& Horiuchill2010l) , weaker AGN jets (e.g. iHondaf 
20091: IrV er et al.l2009h . and galaxy clusters (e.g.l lnoue et al.l2007l : 



Kotera et al 



2009). Candidate transient sources 1 Waxman & Loeb 



2009) include classical Gamma-Ray Bursts ( GRBs; Waxma ni 19951: 
Vietrill 19951 : iMilgrom & Usovlll995l;IWaxmanll2004l2006l : iDermed 
20ld) . low luminosity GRBs (e.g. lMurase et alj2006l) ."AGN flares 
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(e.g. iFarrar & Gruzinovl |2009|). a nd re lativistic ('engine-driven') 
supemovae (e.g. lChakraborti et al 

The arrival directions of UHECRs provide a potentially im- 
portant probe of their origin. Measurements by the Pierre Auger 
Observatory (PAO) rule out isotropy for the hi ghest energy cosmic 
rays at ~ 98% confidence d Armengaudl 120081) . and PAO has fur- 
thermore discovered a correlation between the arrival directions of 
UHECRs with energie s E > 57 EeV and nearby (< 75 Mpc) AGN 
(Abrah am et al. 2008). This result does not, however, imply that 
UHECRs necessarily originate from AGN, because AGN trace lo- 
cal Galactic structure, su ch that the correlation is consistent with a 
variety of other sources dKashti & Waxmarj|2008) ; iGhisellini etail 



120081 : iTakami et all l2009T TTakami & Sato 2009)0 At present the 
sources of UHECR cannot therefore be deduced from their arrival 
directions alone. 

The composition of UHECRs also provides important clues 
to their origin. Although the composition is measured directly at 
low energies (< 10 14 eV), at ultra-high energies it must be inferred 
indirectly by measuring the shower depth at maximum elongation 
Z max . Recent measurements by PAO show that the average shower 
depth (X max ) and its R MS variation decrease systematically moving 
to the highest energies d Abraham et"aiT2 010). This suggests that the 
UHECR composition transitions from being dominated by protons 
below the ankle to being dominated by heavier nuclei with average 
masses similar to Si or Fe at ~ 5 x 10 1 9 eV. We caution, h owever, 
that HiRes has not verified this finding dAbbasi et a l. 2005)0 

The UHECR composition measured by Auger is puzzling. 
One possible explanation is that the accelerated material has an in- 
trinsically 'mixed' composition (with e.g. solar abundances), such 
that protons are accelerated to a maximum energy E = £ p max ~ 
10' 8 ' 5 eV, beyond which only heavier nuclei are accelerated. This 
seems plausible a priori because accelerator size considerations 
show that the maxi mum achievab le energy increases linearly with 
the nuclear charge Z (Hillas 1984). On the other hand, this explana- 
tion appears to require fine tuning because the maximum energy to 
which, for instance, Fe nuclei are accelerated E Fs . m . ix ~ Zx£ p milI ~ 
8 x 10" (Z/26) eV must (by coincidence) be close to the cut-off ob- 
served at ~ 6 x 10 19 eV and expected to occur independently from 
the GZK effect. A 'mixed' composition with metal abundance ra- 
tios similar to the Sun or Galactic cosmic rays also appears incon- 
sisten t with modeling o f the propagation of UHECRs through the 
EBL (jAllard et al. 2008), which suggest that the injected composi- 
tion h as a fairly narrow distribution in charge (e.g. lHooper & Tavlod 

A second possibility is that the accelerated material is dom- 
inated by heavy nuclei. In this case the proton-dominated com- 
position measured near the ankle may be explained as secondary 
parti cles produced by the in teraction of the nuclei with the EBL 
(e.g. lHooper & Tavloj2010l) . A heavy-rich composition is unlikely 
in the case of AGN, galaxy clusters, and supernova shocks because 
the accelerated material originates from the interstellar medium. 
For a solar composition, the fraction of the total mass in Fe nuclei 
and heavier is just X Fc ~ 10~ 3 , such that only for extremely super- 



1 In addition, the latest PAO results suggest that the significance 
of the AGN correlation is reduced from previou s measurements 
iThe Pierre Auger Collaboration: J. Abraham et al. 20q1). 

2 One possible explanation for this discrepency is that PAO and HiRes ob- 
serve different portions of the sky. We note that HiRes a lso does not find the 
correlation of UHECRs with AGN seen by PAO (e.g. Sokols kv & Thomson! 
120071) . 



solar metallicity (~ 10 3 Z o ) could heavy nuclei dominate the total 
UHECR mass. 

In this paper we show that UHECRs from GRBs, unlike AGN, 
may indeed be composed of almost entirely very heavy nuclei. In 
particular, if the outflow from the central engine is strongly mag- 
netized we find that Fe-group nuclei and possibly heavier elements 
(A > 90) are synthesized during its expansion. Although it is well 
established that long du ration GRBs originate fr om the core col- 
lapse of massive stars dWooslev & Blooml l2006t) , it remains de- 
bated whether the central engine is a hyper-accreting black hole 
dWooslevll99lh or a rapidly s pinning, stron gly magnetized neutron 
star (a 'proto-magnetar' ; e.g. lUsov Il992h . We focus here on the 



proto-magnetar model, which recent work has shown can explain 
many of the observed properties of GRBs ( Thompson et al JuOO' 
iMetzger et"al1l2007l; iBucciantini et alj|2007l ; iMetzger et alj|201 
However, similar considerations may apply to accretion-powered 
models, provided that the jet is magnetically-dominated rather than 
a thermally-driven fireball (j j2,2| l. 



2 NUCLEOSYNTHESIS IN MAGNETICALLY-DRIVEN 
GRB OUTFLOWS 

The high temperatures T > 1 MeV near the central engine im- 
ply that all nuclei are dissociated into free neutrons and protons. 
Heavier elements form only once lower temperatures and densi- 
ties are reached at larger radii in the outflow. If the outflow forms 
as a fireball dominated by thermal energy (as would occur if the 
jet i s powered by neut rino annihilation along the rotational axis; 
e.g. lEichleretalJll989h . its entropy is necessarily high 5 > 10 
k\, nucleon -1 . Free nuclei recombine into Helium only once the 
deuterium bottleneck is broken. Since this occurs at low densities 
when the entropy is high, few elemen ts heavier than He are formed, 
simil ar to Big Bang nucleosynthesis dLemoi ne 2002; Beloborodov 
2003). Pure fireballs are therefore unlikely to produce jets enriched 
in heavy elements. 

The situation is different if the jet is accelerated magneti- 
cally, as occurs from proto-magnetars or magnetized accretion disk 
winds. In this case most of the energy is stored in the magnetic field 
(Poynting flux) at small radii and the flow has a much lower entropy 
S ~ 10 - 300 k b nucleon -1 (see Fig. Q] and eq. UJ below). Under 
these conditions Helium recombination occurs at higher densities, 
such that heavier nuclei can be formed efficiently via e.g. the triple- 
a reaction and subsequent a captures. Below we focus on the nucle- 
osynthesis in proto-magnetar winds because the outflow propertie s 
can be calculated with relative confidence dMetzgeretaT]|2010h ; 
however, in 92.21 we briefly discuss the composition of accretion- 
powered outflows. 



2.1 Proto-Magnetar Winds 

When a massive star runs out of nuclear fuel, its core undergoes 
gravitational collapse. This results in a hot 'proto-neutron' star 
(proto-NS), whic h radiates the energy releas ed during the collapse 
in neutrinos (e.g. Bur rows & Latt imer 1986). As neutrinos escape, 
they heat the material above the proto-NS surface, potentially pow- 
ering a supernova (SN) explosion during the first few h undred 
milliseconds after core bounce (e.g. lBethe & Wilsordll985l) . How- 
ever, regardless of how the star explodes, if the core does not col- 
lapse into a black hole, neutrinos continue to heat the proto-NS 
atmosphere on longer timescales t ~ 1 — 100 s. This drives mass 
from the proto-NS into the expanding cavity behind the outgoing 
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Figure 1. Properties of proto-magnetar winds, calculated for a magnetar 
with surface dipole field strength B& p = 5 X 10 15 G, initial spin period 
Po = 2 ms, mass M ns = \AMq, and magnetic obliquity x — 71 1% based on 
the model of Metzger et al. (2010). Wind properties shown include the total 
spin-down power E (solid line), magnetization (To (dotted line), entropy S 
(eq. (3 ; dashed line), and expansion timescale at the radius of He recom- 
bination T exp (dot-dashed line). Also shown is the estimated time required 
for the jet to break out of the stellar surface (1 ~ 10 s), the timescale prior 
to which nuclei are disintegrated by GRB pho tons (t t _n = 10; eq. (6)), and 
the end of the GRB according to the model of lMetzger et al] (2010) (t ss 60 
s). The unshaded area denotes the time interval during which nuclei synthe- 
sized in the wind both survive photodisintegration and may be accelerated 
to energies > 10 20 eV according to the criteria of the internal shock model 
discussed in j j3.2l 



SN shock, producin g 



what is known as a 'neutrino-heated wind' 



dDuncan et al"]| 19861 : burrows et alJI 19951 : lOian & Woosievlll996l : 
hereafter QW96). 

If the proto-neutron star is strongly magnetized (dipole field 
strength Sji p > 10' 5 G) and rapidly rotating (initial spin period 
Pq ~ 1 - 3 ms), i.e. a 'millisecond proto-magnetar', its neutrino 
wind is accelerated primarily by magneto-centrifugal forces. By 
extracting the rotational energy of the neutron star, proto-magnetar 
outflows achieve the power and speed necessary to produce a GRB. 
In our calculations below, we use the time evolution of the power 
E and m ass-loss rate M of prot o-magnetar winds from the detailed 
model o f lMetzgeretal] bOld) . to which we refer the reader for a 
complete description. From E(t) and M(t), we calculate the wind 
magnetization cr (t), which is defined as the ratio of Poynting flux 
to kinetic energy flux at the light cylinder radius m 50(P/ms) 
km. The magnetization equals the Lorentz factor that the jet ob- 
tains if its magnetic energy (Poynting flux) is fully converted into 
bulk kinetic energy. 

Figure [T] shows E(t) (solid line) and <xo(0 (dotted line), cal- 
culated for a magnetar with B Aip = 5 x 10' 5 G, Pq = 2 ms, 
and magnetic obliquity x = f/2 (the angle between the rotation 
and magnetic dipole axes). During the first few seconds after core 
bounce, the magnetar wind is only mildly relativistic (<x < 1) 
because the neutrino-driven mass loss rate is high. However, as 
the proto-NS cools, the outflow becomes increasingly relativistic 
and magnetically-dominated, such that cr reaches > 10 2 - 10 3 on 
timescales ~ 20-50 s. As we describe in ^3] conditions during this 
epoch are ideal for both producing a GRB and accelerating UHE- 
CRs. At late times (t > 100 s), cr increases even more rapidly be- 
cause the neutrino-driven mass loss drops abruptly as the proto-NS 
becomes transparent to neutrinos. Because ultra high-cro have diffi- 



culty efficiently accelerating and dissipating their energy, this tran- 
sition likely ends the prompt GRB emission (Metzger et al. 2010). 
Furthermore, after this point other processes (e.g. y — B or y — y 
pair prod uction) likely tak e over as the dominant source of mass 
loss (e.g. lThompson||2008l) . such that the wind composition could 
change from baryon-dominated to e~/e + pairs. This would appear 
to make UHECR acceleration le ss likely du ring the ultra high-<xo 
phase at late times (although see lAronsll2003h . 



2.1.1 Heavy Element Nucleosynthesis 

The temperature near the proto-NS surface (r = R m > 10 km) 
is set by the balance between neutrino heating and cooling and is 
typically T(R m ) ~ 1 - 2 MeV on the relevant timescales t ~ 10 - 
100 s for producing a GRB and accelerating UHECRs (see QW96; 
their eqs. [46-47]). Because T oc r _1 in the radiation-dominated 
hydrostatic atmosphere, Helium recombination (T = r rcc ~ 0.5 - 
1 MeV) generally occurs just a few NS radii above the surface. 
Because heavy elements begin to form only after recombination, 
the quantity and distribution of nuclei synthesized depends on the 
entropy 5, electron fraction Y e , and the expan sion timescale T exD of 
the outflow at the reco mbination radius (e.g. lHoffman et alii 19971 : 
iMever & Brownll 19971) . where Y e = n p /(n n + n p ) and n p (n n ) is the 
proton(neutron) density, respectively. 

The entropy in proto-NS winds is determined by the amount 
of neutrino heating that occurs in the 'gain region' just above the 
NS surface (but well below the recombination radius), which for 
'normal' (slowly rotating and/or weakly magnetized) proto-NSs is 
well-approximated by the expression 



S(n = 0) = 280C- 1/b L" 1/ V 



1/6-1/3 



v.McV^lO 



2/3 



M14 kbnucleon 



(1) 



where L 51 = L^/10 51 erg s -1 , e v , Me v = %/MeV, R w = R m /W 
km, Mi A = M ns /1.4M G are the electron antineutrino luminosity, 
mean electron antineutrino energy, radius, and mass of the proto- 
NS (QW96; their eq. [48a]); C cs is a correction to the heating rate 
due to inelastic electron scattering (QW96; their eq. [51a]); we 
have included a ~ 20% entropy enhancem ent due to general rel- 
ativistic gravity (e.g. lCardall & Fu ller 1997); and we have assumed 
the electron neutrinos and antineutrinos have similar luminosities 
and mean energies. We calculat e L v (t), e v (t ) , and R m (t) using the 
proto-NS cool i ng cal culations of lPons et al.l i 19991) as described in 
iMetzgeret ail ( 1201 fj) . 

The entropy in proto-magnetar winds depends on the obliq- 
uity angle x- For aligned rotators (x ~ 0), material leaves the NS 
surface near the rotational pole and its entropy is similar to the non- 
rotating value, i.e. S(Q.,x = 0) ~ S(C1 = 0) (eq. (T|). By contrast, 
oblique rotators (x ~ it 12) lose most of their mass in outflows from 
the rotational equator. Equatorial outflows experience significant 
centrifugal acceleration in the gain region, which reduces the heat- 
ing received by the outflowing material and suppresses the entropy 
exponentially, viz. 



S(Q.,x = x/2) = S(£l = 0) x exp[-P ccM /P], 



(2) 



where P c , nt ss 2. IL. 



■„ ms is taken from the numerical calculations 
of lMetzger et alj (2007) for a proto-NS with M ns = 1.4M Q and ra- 
dius R m = 10 km (their eq. [37]). Although the dependence of P C!:ni 
on R ns has not yet been determined, the proto-NS has already con- 
tracted to its final radius =s 10 km by the times of int erest for accel- 
erating UHECRs (t > 10 s), such that the results of Metzg eret al.l 
( 120071) are applicable. 

The expansion timescale in the wind is defined as r cxp = 
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('"/v r )lr rcc , where v r is the outflow velocity and we calculate the 
recombination temperature T Iec from the density and entropy, as- 
suming nuclear statistical equilibrium (e.g. QW96; their eq. [62]). 
We calculate the velocity using mass continuity M = pv r A, where 
A = 4-7iR^ s fopea(r/R ns ) 3 is the areal function of the dipolar flux tube 
(recombination generally occurs interior to the light cylinder) and 
/open ~ ^ns/2^L is the fraction of the magnetosphere at the surface 
open to outflows. 

The entropy S(t) and expansion timescale T,, xp (f) for the wind 
solution in Figure [T] are shown with dashed and dot-dashed lines, 
respectively. Note that S is approximately constant in time because 



the rising value of S(Q, = 0) oc L- 



(eq. (T|) is offset 



by an increase in the exponential factor P ccn JP (eq. 0). Also note 
that at late times (> 100 s) the expansion timescale T cxp becomes 
comparable to the timescale over which the properties of the wind 
are changing, suggesting that the steady-state assumption we have 
adopted may break down. Although this does not affect the con- 
clusions of this paper because we are focused on the jet properties 
on timescales of tens of seconds, future numerical work is required 
to more accurately address the properties of magnetized proto-NS 
winds at late times and low neutrino luminosities. 

The electron fraction Y e is important for two reasons. First, 
the value of Y e determines the channel by which helium burns to 
form carbon. Since this is the slowest reaction, Y e impacts the total 
heavy element yield. Under proton-rich conditions (Y e > 0.5), 12 C 
forms via the standard triple-a reaction sequence 4 He(2a',y) 12 C. 
For Y e < 0.5, on the other hand, the n eutron capture chan- 
nel 4 H e(ttn,y) 9 Be(«,n) l2 C instead dominates iWooslev & Hoffmanl 
1992). The electron fraction also determines the distribution of 
heavy nuclei synthesized. When Y e < 0.5, nuclei up to the 
N=50 neutron closed shell (A ~ 90) are created via alpha 
particle and ne utron captures, depending on Y e and the final 



a fraction (e.g. IWoosley & Hoffmanl 1 19921 : iRoberts et all l20ld : 
lArcones & Montesl l2010h . If Y, is sufficiently low, even heav- 
ier r-process elements (with characteristic peaks at A w 130 
and A =s 195) can b e created by additional neutron captures 
(e.g. lSeegeret"aT1ll965h . Under proton-rich conditions (Y e > 0.5), 
by contrast, mainly Fe-group elements (A ~ 40-60) are created, al- 
though for Y e > 0.55 elements up to A « 64 may be created by pro- 
ton captures before the flow reaches 'bo ttleneck' nuclei with long 
/3-decay timescale s such as 64 Ge (e.g. lArcones & Monies 2010; 
IRoberts et aifcoid) . 

Although the nucleosynthesis is sensitive to Y e , its value in 
proto-NS winds is rather uncertain. The surface of the proto- 
magnetar is neutron rich (y c < 0.1), but as nucleons are accelerated 
outwards in the wind they are irradiated by electron neutrinos and 
antineutrinos. This drives Y e to a value ~ 0.4 - 0.6 that depen ds 
on the precise v e /v e luminosities and spectra JOian et al . 1993)0 
The value of Y e is thus sensitive to the details of neutrino trans- 
port and interactions (e.g. lRampp & Jan ka 2000; Mezzacappa et al.l 



3 Neutrino absorptions may allow these bottlenecks to be circumvented un - 
der some circumstances (the so-called 'v— p process ': lFrohlich etai]|20O6l) . 
However, this is unlikely to be relevant in proto-magnetar winds because the 
e xpansion timescale is much shorter than the neutrino capture timescale. 

4 iMetzger. Thompson & Ouata ert ( 2008) found that Y e can remain low if 
the proto-magnetar is rotating extremely rapidly (P < 1 ms), such that mat- 
ter is accelerated away from the surface before its composition is altered by 
neutrinos. However, at the late times of interest for accelerating UHECRs 
the proto-NS is unlikely to be rotating this rapidly because it loses appre- 
ciable angular momentum to winds during even the first few seconds after 
forming. 



1200 ll ; iDuan et alj|201(il) . Early SN calculations found that Y e de- 
creased from > 0.5 a t early times to < 0.5 as the proto-NS cooled 
JWooslevetliflll994l) . More recent calculations including addi- 
tional neutral-current interactions, however, find that Y e rises from 
~ 0.5 to ~ 0.6 on a timescale ~ 10 seconds (Hudepohl et alj|2010h . 
Unfortunately, few calculation s have yet been perform ed in the case 
of rapidly spinning proto-NSs. Thomp son et al.l ((2005) find that the 
v e /% luminosities and temperatures at t < 0.6 s are appreciably al- 
tered in the case of core collapse with rapid rotation, such that Y e 
is driven to a higher value than in the non-rotating case. It is, how- 
ever, difficult to extrapolate their results to the much later times 
~ 10 - 100 s of interest here. 

In our calculations below we consider both possibilities Y e < 
0.5 and Y e > 0.5, keeping in mind that the true electron fraction 
probably lies in the range 0.4 < Y e < 0.6 and could vary with 
time and from event to event, depending on e.g. the NS mass, rota- 
tion rate, and obliquity. However, given the uncertainties, we can- 
not rule out the possibility that Y e is appreciably lower, in which 
case even heavier r- process elements with A > 100 are produced. 
In particular, although present theory does not find the neces- 
sary conditions for the second or third peak r-process in proto-NS 
winds, indirect evidence (from e.g. Galactic chemical evolution) 
suggests that the r-process in deed originates from core-collapse 
SNe (e.g. lMathews et al.lll992l) . 

Given 5, T cxp , and Y e , we estimate the total mass fraction X n 
synthesized in h eavy nuclei (A > 56) using the following analytic 
expressions from Roberts et al. I d20 1 Oh (their eqs. [B3] and eq. [B3] 
and rem. cf. lHoffman et al.lll997l) : 



{l-exp[-8xl0^ 3 (^)( E ^ r )- 3 ]}, 



(3) 



_ 2l -l/2 



Y e < 0.5 
Y. > 0.5 



(4) 



The critical quantities in square brackets are essentially the reac- 
tion rate for 4 He — > 12 C integrated over the timescale available 
for burning At ~ r exp . In the proton-rich case, the reaction rate 
depends on the entropy squared because triple-a is an effective 
three-body reaction and p oc 1 /S at fixed temperature, while in the 
neutron-rich case the reaction rate depends on entropy cubed be- 
cause 4 He(cm,y) 9 Be(a',n) 12 C is an effective four-body reaction. 



2.1.2 Results 

Figure|2]shows the total mass fraction in heavy nuclei X n (eq. |4]) 
as a function of time after core bounce, calculated for a magnetar 
with B dip = 5 x 10 15 G and P = 2 ms, as in the wind model shown 
in Figure Q] We show four models, corresponding to different val- 
ues of the magnetic obliquity ;^ and the wind electron fraction Y e . 
In the case of a neutron-rich outflow (Y e = 0.4), both aligned and 
oblique rotators have X n ~ 1 throughout the entire wind phase; 
heavy elements form efficient ly because the 4 He(«n,y) 9 B e(o',n) l2 C 
channel is fast and available dWooslev & Hoffmanll 19921) . By con- 
trast, proton-rich outflows have < 1 at all times because carbon 
must form through the (slower) triple-a channel. Note that at fixed 
Y e = 0.6, the heavy fraction is larger for an oblique rotator because 
X], depends sensitively on the wind entropy, which is suppressed in 
equatorial outflows (eq. |2"|). Although the precise value of X h de- 
pends on the (uncertain) value of Y e , we conclude that in all cases 
a significant fraction of the mass of the wind is locked into heavy 
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Figure 2. Fraction of the mass in proto-magnetar winds that is synthesized 
into heavy elements X b (eq. fill, calculated for a magnetar with dipole field 
strength B L \ lp = 5 X 10 15 G, initial spin period Pq = 2 ms, mass M ns = 
IAMq (see Fig. [TJ. Four models are shown, calculated for both aligned 
(magnetic obliquity x = 0) and oblique (x = ff/2) rotators and for neutron- 
rich (electron fraction Y e = 0.40) and proton-rich (Y e = 0.60) outflows. The 
area shown in white denotes the epoch during which nuclei synthesized in 
the wind may be accelerated to energies > 10 20 eV according to the internal 
shock model discussed in j )3.2l 



nuclei during the epoch of UHECR acceleration. Whatever mass is 
not used to form heavy nuclei remains as Helium, i.e. Xn s — 1—Xh- 
Although heavy nuclei may form deep in the wind, they 
could in principle be destroyed at larger radii in the jet. In par- 
ticular, spallation can occur if heavy nuclei collide with a parti- 
cle of relative energy exceeding the nuclear binding energy as 8 
MeV nucleon -1 . For heavy nuclei, the cross section for inelas- 
tic nuclear collisions is ~ 1 barn, similar to the Thomson (pho- 
ton) cross section. Because the flow contains ~ Z electrons per 
nucleon, this implies that inelastic scattering is only important 
at ra dii well inside the Tho mson photosphere at r p h ~ 10 11 cm 
(e.g. iKoers & Gianniosll2007l) . If strong shocks were to take place 
at these small radii, the resultant heating could in principle destroy 
the nuclei. However, for the strongly magnetized jets under consid- 
eration, the magnetization is probably still high at these radii, such 
that s trong shocks are highly suppressed (e.g. iKennel & Coronitil 
1 1984) . In particular, the radial profile of acceleration to a Lorentz 
factor T oc r 6 occurs more gradually in MHD jets (0 



1/3 - l/2;|Prenkhahn & Spruitl2002l:lviahakis & Konigill2003allbl: 
Komissarov et alj 120091 ; iTchekhovskov et alj 120091 ; iGranot et al.l 
20ld) than in thermally-accelerated fireballs (j8 = 1). This implies 
that full acceleration to T ~ o"o ~ 100 - 1000 is generally only 
possible at large distances (typically > 10 12 cm), well into the re- 
gion where the plasma is collisionless with respect to direct nuclear 
collisions. We thus conclude it is unlikely that nuclei are destroyed 
during the collimation and acceleration phase of the jet (although in 
93.3| we discuss the conditions under which nuclei are destroyed by 
GRB photons during the subsequent UHECR acceleration phase). 



2.2 Accretion-Powered Outflows 

Heavy nuclei may also be produced in GRB outflows powered by 
black hole accretion, provided that the jet is magnetically-driven. 
If the jet threads the black hole event hori zon (as occurs if its 
powe r derives from the spin of the black hole; Blandford & Znajek 
1 19771) . then the outflow composition is effectively baryon-free near 



the ergosphere. To what degree baryons are entrained in the jet 
at larger radii depends on uncertain diffusive processes from the 
jet walls, thereby making cr difficult to predict with confidence 
(e.g. iLevinson & Eichlerll2003l: lMcKinney||2005l) . The entropy in 
this case depends on the amount of heating (due to e.g. v—v annihi- 
lation or magnetic reconnection) and this uncertain baryon loading. 

If, on the other hand, the jet directly threads the sur- 
face of the accretion disk (e.g. Blandford & Pavnel Il982l) . its 
mass-loading may (as in proto-magnetar winds) be controlled 
by neutrino heating in the disk atmosphere (e.g. iLevinsonl 
20061: ISurman et aljhood iMetzger. Thompson & Ouataertl 1200^ 



Metzg er. Piro & O uataert 2008). In this case, the outflow i s likely 



to be proton-rich IMetzger. Thompson & Ouataertl [2008). with 
an entropy and expansion timescale similar to those in proto- 
magnetar winds, but depending on the accretion rate, black hole 
mass, o pen magnetic field geometry, an d wind launching radius 
(see e.g. IMetzger. Piro & Ouataert 2008 ; their eqs. [24] and [26]). 
Although the outflow properties are less certain than in proto- 
magnetar winds, efficient nucleosynthesis may well lead to a heavy 
fraction ~ 1 in accretion-powered outflows as well. 



3 UHECR ACCELERATION IN PROTO-MAGNETAR 
JETS 

Near the light cylinder radius ~ 100 km the power in 
the magnetar wind is con centrated in the rotational equator 
(e.g. iBucciantini et alj l2006h . On larger scales the wind is col- 
limated i nto a bipolar jet by its interac t ion with the progenitor 
star (e.g. lUzdenskv & MacFadvenl 120071 ; IBucciantini et alJl2007t 
120081 12009). After the jet propagates through the star and 'breaks 
out' of the surface on a timescale t ~ 10 s (e.g. lAlov et all 
2000), the magnetar wind escapes through a relatively clear chan- 
nel. High energy emission (the 'GRB') occurs when the jet dissi- 
pates its energy via shocks or magnetic reconnection at larger radii 
~ 10 13 - 10 16 cm (see below). Although there are many potential 
sources of short timescale variability in the outflow (e.g. interac- 
tion of the jet with the confining stellar envelope), numerical sim- 
ulations show that the time- and angle-averaged values of E(t) and 
o"o(f) (Fig.Q3 mat ch those set by the magnetar wind at much smaller 
radii (e.g. Morsony et al. 2010). We assume that the opening angle 
of the jet at the stellar surface is 6 S ~ 4° (e.g. IBucciantini et al.l 
2009), such that the isotropic jet luminosity is related to the wind 
power by E lso = Ej f b , where f b = dj/2 = 2x 10~ 3 is the beaming 
fraction. 

If UHECRs originate from GRBs, they are probably acceler- 
ated by the same dissipative mechanisms responsible for accelerat- 
ing electrons and powering the GRB emission. For magnetically- 
dominated outflows, most of the jet's energy resides in Poynting 
flux near the central engine; a sizable fraction of this magnetic en- 
ergy must be converted into bulk kinetic energy in order to explain 
the hi gh Lorentz factors (r > 10 2 ) inferred from GRB observations 
(e.g. lLithwick & SarfcOOll) . Depending on the means and efficacy 
of acceleration in the jet, gamma-ray emission (and UHECR ac- 
celeration) may be powered either by the dissipation of the jet's 
Poynting flux directly ('magnetic reconnection') near or above the 
photosphere ( 93. lh and/or via 'internal shocks' within the jet at 
larger radii Q3.2t . 

Below we discuss the conditions for UHECR acceleration in 
magnetar jets. Our treatment closely follows past work that relie s 
on direct constraints from GRB observations (e.g. I Waxmanll 1 995T) . 
However, this is the first time that the analysis has been applied 



© 2010 RAS, MNRAS 000,[T|-?? 



6 B. D. Metzger, D. Giannios, & S. Horiuchi 



Magnetic Reconnect 'on 




10"' 



20 30 40 50 60 

Time (s) 

(a) Magnetic Reconnection 
Internol Shoc<s 




20 30 
Time (s) 

(b) Internal Shocks 

Figure 3. The maximum energy £ lm , to which Fe nuclei can be accelerated 
in the jet (left axis) as a function of time after core bounce, calculated for the 
wind solution from Figure[T]and shown for both reconnection (top panel; 
jj3.1> and internal shock (bottom panel; j|3.2) models. In each case both 



the dynamical timescale constraint ? ac 



and the synchrotron cooling 



constraint ? acc < ? coo i are shown with solid lines; the limiting (minimum) 
value of £" max at each time is shown in boldface. Dotted lines show the 
optical depth to nuclear disintegration via the GRB photons r y -fj (right 
axis; eq. (6j). The unshaded area denoted 'HEAVY UHECRs' represents 
the epoch during which nuclei are simultaneously capable of reaching ultra- 
high energies (i.e. E mdx > 10 20 eV) yet are not destroyed by GRB photons 
(i.e. T y _,v < 10; see eq. (6]). 



using a 'first principles' central engine model for the jet properties 
and the location where dissipation occurs. 



3.1 Acceleration by Magnetic Reconnection 

One way reconnection can occur in the jet is if the rotation and 
magnetic axes of the NS are misaligned (x > 0), such that the out- 
flow develops an alternating or 'striped' magnetic field g eometry 
on the scale of the light cylinder radius R L dCoronitilll990h . A sim- 
ilar field geomet ry may result if the jet is susceptible to magnetic 
instabilities (e.g. [Giannios & Sp ruit 2006; Mc Kinney & Blan dford 
1 20091 ; Moll 2009). If this non-axisymmetric pattern is preserved 
when the flow is redirected along the polar jet, the resulting geom- 
etry is conducive to magnetic reconnection. We adopt the model de- 



veloped bv lDrenkhahn & Spruiu J2002h . in which magnetic dissipa- 
tion occurs gradually from small radii up to the 'saturation' radius 
>" ~ Rrn-ig — o-qPc/6£, beyond which reconnection is complete and 
the flow achieves its terminal Loren tz factor, where e = vjc and 
v r is the reconnection speed (see also Lyubarsky 2010). During this 
process, approximately half the Poynting flux is directly converted 
into kinetic energy (producing acceleration) and the other half is 
deposited into the internal (thermal) energy of the flow. We assume 
a fixe d value e = 0.01 independent of radius (e.g. lUzdenskv et all 
2010), but our results would be qualitatively unchanged if recon- 
nection were 'triggered' abruptly when, f or instance, the outflow 
transi tions to a collisionless regime (e.g. iMcKinnev & Uzdenskvl 
l20ld) . 

One important constraint on a potential UHECR source is 
the maximum energy £ raax to which cosmic rays can be acceler- 
ated. One mechanism for accelerating particles in regions of mag- 
netic reconnection is first-order Fermi acceleration, which occurs 
when particles a re deflected at th e converg ing upstream in the re- 
connection flow (Giannios 2010; see e.g. Lvutik ov & Ouvedll2007l 
for an alternative possibility). For sufficiently fast reconnection, 
the acceleration timescale is very short, similar to the Larmor gy- 
ration timescale f acc ^ t]. X{: 2kE' /ZeB'c, where E' = E/F and 
B' m (Z? iso r~ 2 r~ 2 c~') 1/2 are the particle energy and magnetic field 
strength (evaluated at r = R mal! ), respectively, in the jet rest frame; 
T ~ <To/2 is the bulk Lorentz factor in the acceleration zone; and 
'/acc ~ 1 is a fudge factor that parametrizes our ignorance of e.g. the 
precise reconnection geometry. One constraint on £ max is that f acc 
must be < the jet expansion timescale ? exp =; R mag /Fc, such that cos- 
mic r ays can be ac celerated within the dynamical timescale of the 
flow dHillasH 19841) . If the expansion timescale constraint does in- 
deed determine £ max , heavy nuclei achieve a value of £ max that is a 
factor of ~ Z larger than protons. Another constraint arises because 
the acceleration competes with cooling of the nucleus. The domi- 
nant cooling mechanism is synchrotron emission, which occurs on 
a timescale 



3A 4 mp c 7 r 
E 



(l0 20 ev) 



g (se) (2s) > 
100) I 



10 52 erg s 



10 13 cm/ 



(5) 



where e mag is the fraction of jet power carried by Poynting flux and 
is ~ 1 in the case of reconnection-powered outflows. 

The top panel of Figure [3] shows £ max in the reconnection 
model as a function of time, based on the two independent con- 
straints f acc < f cxp and f acc < t cool and calculated for the same mag- 
netar wind shown in FigureQ] assuming rj X{: = 1. We also assume 
a pure Fe composition, although our results are similar if the com- 
position is instead dominated by heavier nuclei (A > 90). Note that 
at the earliest times (f < 15 s in this example) synchrotron losses 
place the most severe constraint on E miix , but at later times (t > 15 s) 
the expansion timescale constraint is more severe. As we discuss in 
93.3| below, a more severe constraint at early times arises because 
heavy nuclei can be disintegrated by the GRB photons. The time 
interval shown in white (20 s < t < 50 s) denotes the epoch during 
which heavy nuclei both survive photodisintegration and achieve 
values of £ max > 10 20 eV which are sufficient to explain UHECRs. 
During this epoch ~ 10 51 ergs of rotational energy is extracted from 
the magnetar, a large fraction of which could be placed into UHE- 
CRs. 
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3.2 Acceleration by Internal Shocks 

If the jet is accelerated efficiently, then a significant fraction of its 
Poynting flux is converted into kinetic energy. Because cr ~ T in- 
creases monotonically during the GRB (Fig. QJ, slower material is 
released prior to faster material. Strong shocks will occur once the 
faster material catches up, provided that the r esidual magnetiza- 
tion o f the jet at the collision radius is < 0.1 dKennel & Coronitil 
1 1984b . This scenario is similar to the standard internal shock 
model for GRB emission (e.g. Ree s & Meszarosl 1 19941) . except 
that the shock occurs between the fast-moving jet and the 'bulk' 
shell composed of the total mass re leas ed at earlier times ( see 
Metz ger. Ouataert & Thomp son 2008 and Metzger et al. 2oTol for 
a detailed description). 

If UHE CRs are accelerated by internal shocks (e.g. via Fermi 
acceleration; iGallant et alJI 19991) . the conditions on £ max are quite 
similar to the case of magnetic reconnection discussed above 
(i.e. fa,. c < f cxp and f acc < f cooI , with 7/ acc ~ 1), except that (1) shocks 
occur at larger radii; (2) the Lorentz factor of the shock is no longer 
the instantaneous value of the jet T ~ erg, but rather the average 
of the material released since jet break-out; and (3) the magnetic 
field necessarily carries a smaller fraction e B «: 1 of the jet power 
than in the reconnection case (e B ~ 1), such that strong shocks are 
possible. In the bottom panel in Figure [3] we show the constraints 
on £ max m the internal shoc k model, assuming e B = 0.1, 7/ acc = 1 
(e.g. lAchterberg e t al. 2001), and an Fe composition. Note that syn- 
chrotron losses are unimportant in this case, such that the expansion 
timescale constraint t. xc < r cxp sets E max at all times. Again, the time 
interval in white denotes the epoch when heavy nuclei survive pho- 
todisintegration (see below) and £ max > 10 20 eV and ending only 
when the GRB ends at / w 60 s. We conclude that regardless of 
whether jet dissipates its energy via reconnection or shocks, UHE- 
CRs with £ max > 10 20 eV can be produced simultaneous with the 
epoch of GRB emission. 

3.3 Photodisintegration at the Source 

For heavy nuclei to escape the jet, they must avoid photodisinte- 
gration (primarily via giant dipole resonances; GDR) and other en- 
ergy loss mechanisms (e.g. pion production, two nucleon emission, 
photo-absorption) due to interaction with the GRB photons. The 
'Band' function that characterizes GRB spectra peaks at a char- 
acteristic energy E p — 0. 1 — 1 MeV. This corresponds to an en- 
ergy E p as E p /T in the rest frame of the jet and an energy E p « 
E p Ja/T 2 ~ 60 GeV^p/SOOkeV^oA^r/lO 2 )- 2 in the rest frame 
of a cosmic ray with (observer frame) energy E = y^AmpC 2 , where 
£20 = .E/10 20 eV and A 56 = A/56. GDR and pion production occur 
at typical energies ?gdr ~ 18Aj° 21 MeV and e A ~ 0.3 GeV, respec- 
tively, which are both typically below the Band peak E p . Below the 
peak the Band spectr um has a power-la w shape dN y /dE y oc E y , 
where a a -1 (e.g. Na va et ail l2010h . For a = -1 the num- 
ber of photons per logarithmic frequency interval around the res- 
onance is approximately constant, such that in the jet rest frame 
N y = j(dN y /dE y )dE y * CU y /E' p , where U y « £ iso e rad /(47rr 2 C r 2 ) 
is the total photon energy density, e rad ~ 0.1 - 1 is the radiative 
efficiency of the jet, and C ~ 0.2 is the fraction of the gamma-ray 
energy released below the Band peak. 

From these considerations we estimate that the number of in- 
teractions (photodisintegrations or pion-producing) experienced by 
an UHECR is given by 

r y _ w » N 7 cr t (Ae r / et )AR * — , (6) 



where AR is the pathlength traversed by a nucleus as it is accel- 
erated and escapes the acceleration region, which is approximately 
equal to the characteristic radius of the outflow in the comoving 
frame x r/T; and <r r , e r and Ae r are the resonance cross-section, 
energy, and width, respectively. For heavy nuclei GDR dominates 
pion production and other energy loss mechanisms, so we assume a 
line width Ae, /e r ~ 0.4A!? 6 21 and cross section cr r ~ 8xlQ~ 26 A5 6 cm 2 
appropriate for G DR (e.g. iKhan et alj 120051 : iMurase et all 120081: 
iPe'er et alj2009l) . 

On the right axes in Figure [3] we plot r y - N (eq. |6]) as a func- 
tion of time in both reconnection and internal shock models, cal- 
culated assuming e rad = 0.5, C = 0.2, and a pure Fe composition. 
We define the criterion for the survival of a nucleus as r y - N < 10 
because a nucleus can experience ~ 10 photodisintegrations before 
its composition is appreciably altered (see also 34. 11 below). Note 
that T r _iv 3> 10 immediately after jet break-out, which shows that 
heavy nuclei synthesized at early times are destroyed. However, be- 
cause T and the acceleration radius r both increase monotonically 
during the burst r y - N oc r~ 2 r~' decreases rapidly with time, such 
that heavy nuclei survive (i.e. r y - N < 10) at times t > 15 - 20 s. 



4 FROM THE SOURCE TO EARTH 

The propagation of UHECRs through intergalactic space is impor- 
tant because their attenuation by the CMB/EBL and deflection by 
magnetic fields determines the maximum distance of sources and 
the time delay tj between the GRB and their arrival on Earth. Al- 
though the propagation of proton s and nuclei have been studied 
extensively in previous works (e.g. | AUard et al.l20 05 ; Taka mTet al.l 
120061: lAllard et al.ll2008l ; lHooper & Tavlor!l2010l) . we extend these 
considerations to the novel case of ultra-heavy nuclei. 

4.1 Photodisintegration by the EBL 

In this section we estimate the distance a heavy nucleus travels 
before it loses a significant fraction of its energy to interactions 
with the CMB and EBL (neglecting the effects of magnetic fields 
for the moment; in 94.21 we discuss the validity of this assump- 
tion). As in the GRB jet, Giant Dipole Resonances (GDRs) are 
the most important interaction. A nucleus of mass A and energy 
Ek = yAAm p c 2 interacts via GDR with a EBL photon of en- 
ergy 6q DR ^ Egdr/Ta - 5 x 1Q~ 3 A 56 E~q eV. Nuclei with energies 
E ~ io 19 - 5 ~ 20 eV and masses less than that of Fe generally inter- 
act with the CMB, which dominates over other photon sources at 
E < E C mb ~ 5 x 10~ 3 eV. Intermediate mass elements (e.g. He, 
C, O, Si) with ultra-high energies therefore have very short path- 
lengths through the EBL and are like l y to arrive at Eart h as sec- 
ondary protons (e.g. lAllard et al.ll2005l;lAllard et alj2008h . 

Nuclei in the Fe group or heavier with E 2 q ~ 1, by con- 
trast, have e GDR > £cmb and instead interact with the the extra- 
galactic infrared background ( EIB). We approximate the EIB spec- 
trum as dN y /dE y oc E~ 2 - 5 (e.g. lMalkan & SteckeJll998l;lDole et all 
2006), such that the number of photons per decade in frequency 
N y = J(dN y /dE y )dE y oc Ez 15 . Thus, for nuclei with £q DR > E C mb 
(i.e. A^E^q > 1), we (crudely) estimate the mean free path for a 
single GDR interaction as 

-lo - — -, - 1 — — « WE£ 5 A° s I Mpc, (7) 

N y (E = e GDR )cr GDR (Ae GDR /e GDR ) 

where we have used the fact that c tgdr k A, 6qdr <x A ~ 2 , and have 
normalized Aq using the results of lAllard et alj fc005t) . 
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Equation l[7} represents the distance before a single nucleon 
or a particle is ejected. In order to calculate the total energy loss 
length xi5 (defined as the distance until a nucleus loses ~ 25% of 
its initial energy), we must determine how A changes as the nucleus 
is torn apart. Because nuclei are ejected at subrelativistic speeds in 
the rest frame of the parent nucleus, in the lab frame they share 
the same Lorentz factor as the parent. Thus, the parent energy E 
decreases proportional to its mass oc A, with its Lorentz factor re- 
maining constant at its initial value y A . As a result, the parent nu- 
cleus continues to interact with the same radiation background of 
energy ~ £gdr/7a, i-e. N y ~ constant. However, because the disin- 
tegration cross section decreases (A6 GDR /6 G dr) X o"gdr k A 1 2 , the 
mean free path increases as A oc A~ L2 . The total energy loss length 
is therefore given by the sum of the individual mean free paths 

A A „A 

Xl5 * n;/ J] A - = "ej'-loA 1 " 2 £ r u « A A i2 h^ x- [2 dx 

i=3A/4 i=3A/4 J3A/4 

~ 03A Q Ah- c] ' » nOE^ 3 A^n c] -' Mpc, (8) 

where ~ 1 is the mean number of nucleons ejected per disinte- 
gration. 

Equation ([8} illustrates that at a fixed (measured) energy E ~ 
10 20 eV, the distance of accessible sources increases as D ~ xis K 
A . Although the accessible distance to Fe-rich sources is (by co- 
incidence) similar to that of protons (xis ~ 150 Mpc), sources rich 
in nuclei with A « 90(200) nuclei (as predicted if the jet is neutron- 
rich) are observable to a distance ~ 2(5) times larger. Such an in- 
crease in the number of accessible volume potentially alleviates en- 
ergetic constraints on GRBs as UHECR sources (e.g. lEichler et al.l 
120101) by expanding the number of sources within the GZK horizon. 
We emphasize, however, that our above analysis has neglected sev- 
eral potentially relevant details (e.g. evolution of the EBL with red- 
shift, energy losses to e~/e + production). A full propagation study 
will be necessary to determine whether a composition rich in ultra- 
heavy nuclei X h ~ 1 (and with a potentially significant He mass 
fraction Xn c — 1 - Xh) is consistent with the measured UHECR 
spectrum and composition. 

4.2 UHECR energetics and limits on the intergalactic 
magnetic field 

We now discuss constraints on GRB energetics, and place limits on 
the strength of the intergalactic magnetic field, that are consistent 
with the hypothesis that heavy nuclei from GRBs are the primary 
source of UHECRs. 

If the injected UHECR spectrum is a power-law dN/dE oc 
E~ s with s > 2, then the minimum source power per volume 
£ m i„ required to explain the observed UHECR flux is dominated 
by the smallest observed energies ~ 10 19 eV and is given by 
Smm ~ E 2 § |, nl9 v « (0.5-2) x 10 44 erg Mpc V' 1 depend- 
ing on j dWaxmanlll995l ; lBerezinskv et 31.1 2006; Kat z et alj|2009l ; 
M urase & Ta kami 2009), although this value may be somewhat 
smaller for a composition containing ultra-heavy nuclei because 
of their larger accessible distance oc xis (eq. fU). Each GRB 
must then on average contribute an (isotropic) energy 6 C r£q^ b > 
10 54 (//10)(JV G rb/0.5 Gpc~ 3 yrV ergs, where e CR is the ratio 
of the energy placed into UHECRs to that radiated as ga mma- 
rays; JV GRB is the local GRB rate (e.g. iGuetta et all l2005h : and 
/ = f E(dN/dE)dE/& min > 1 is the factor by which £ min under- 
stimates the total cosmic ray energy. For s > 2, the correction / 
is generally large (> 10), which may place severe constraints on 
GRB models for UHECR by requiring that ecr be unphysically 



high (e.g. lEichler et alJlzoTob . If, on the other hand, the acceler- 
ated spectrum is shallow (s < 2), the total energy is dominated by 
the largest energies ~ £ max , such that the lower limit on ^cr^gm' ls 
probably much less severe. As we discuss in ij5] a shallow injected 
spectrum also appears to be consistent with several other features 
of GRB models dominated by heavy nuclei. 

We now discuss constraints on the intergalactic magnetic field. 
Although the (undeflected) path length of heavy UHECRs may be 
quite large (eq. (8)), heavy nuclei can be deflected significantly by 
the intergalactic field. In particular, too strong of a field results in a 
large deflection angle Q d > 1 rad, which would signi ficantly reduce 
the accessible volume of sources (e.g. IPirarfo lO). For an inter- 
galactic magnetic fields B- mt of correlation length 6, the deflection 
angle of a particle traveling a distance D oc X75 is 

0i ~ ^(»/<5) 1/2 <5/fl,, (9) 

where R t = E/ZeB iut is the Larmor radius of the particle. Setting 
Q c i ~ 1 rad gives an upper limit on the magnetic field strength: 

Aim < 7 x lO-^^oD^^Xo'G, (10) 

where A 00 = D/100 Mpc, S Mp c = 6/1 Mpc, and Z 30 = Z/30. 
This limit is not particularly stringent because observational probes 
of intercluster magnetic fields constrain their strength on cor- 
relation lengths 6 ^ , 1 M pc to be less than ~ InG (see, e.g., 
iNeronov & Semikozj|2009h . 

A minimum intergalactic field strength can be derived by the 
requirement that ~ 10 20 eV cosmic rays suffer the minimum de- 
flection that results in a spread in their arrival time after the GRB, 
such that at least one source within the GZ K volume contributes t o 
the local UHECR flux at any instant (e.g. lWaxman & Loebll2009l). 
For a n observed local rate A/qrb ~ 0.5 Gpc~ 3 yr~' ( IGuetta et al.l 
120051) . the rate that GRBs occur within a sphere of radius D is 
H ~ 2 x 10~ 3 yr~'£> 3 00 . Since the duration of a GRB ? G rb is « 
the cosmic ray signal must reach Earth with a spread in arrival times 
td > K~ = 500Z)7Q yr, such that at least one 10 20 eV source con- 
tributes at any given time. Since the time delay is related to the 
deflection angle via tj ~ D6^/4c, we use equation ® to arrive at a 
lower limit on the magnetic field strength: 

S int > 2 x lQr u E7DD$?6£gz£Q, (11) 

For A ~ 100 and Z ~ 40 nuclei the attenuation length is D ~ xis ~ 
300 Mpc (eq. QQ), resulting in the requirement B mt > 10 -13 G. 

Although the strength of the extragalactic magnetic fields is 
poorly constrained by observations, the limits given by equations 
d 10| > and 1 11 It allow for a reasonably wide range of field strengths 
consistent with heavy nuclei from GRBs as the source of UHE- 
CRs. We note that, in addition to the intergalactic field, UHECRs 
may interact with the more localized fields associated with large 
scale structure formation (e.g. galaxy clusters). Although it is dif- 
ficult to predict the magnitude of this effect with confidence, the 
deflection angle and delay time could in principl e be significantl y 
enhanced (e.g. iDas et alTl2008l : see, however, iDolag etalj|2005h . 
Likewise, a constraint similar to equation 1 11 11 1 arises by requir- 
ing that the angular distribution of sources on the sky be consis- 
tent with with the lack of small scale anisotropy observed by PAO 
(Taka mi & Satol l2008); evaluating this constraint in the context of 
ultra-heavy nu clei is, however, nontriv ial and will require addi- 
tional work (cf. lMurase &~T akami 2009h . 
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5 DISCUSSION 

The spectrum and composition of UHECRs measured by Auger 
are consistent with an accelerated composition dominated by heavy 
nuclei. However, until now there has been no astrophysical motiva- 
tion for considering such a source. In this paper, we have shown 
that magnetically-dominated GRB outflows may synthesize heavy 
A a 40 - 200 nuclei as they expand away from the central engine 
(Fig. 0. Focusing on the millisecond proto-magnetar model, we 
have shown that the regions of magnetic reconnection or shocks 
responsible for the GRB emission also allow heavy nuclei to be 
accelerated to ultra-high energies £ raax > 10 20 eV while not being 
disintegrated by GRB photons (Fig. [3}. 

Models that invoke GRBs as the source of UHECRs have 
been criticized on the grounds that the required energetics may 
be insufficient if the rate and observed gamm a-ray fluences from 
GRB s are a proxy for th e UHECR flux (e.g. iFarrar & Gruzinovl 
120091 : lEichler et al1l2010l) . However, this argument depe nds on the 
(uncertain) local rate of GRBs (e.g. iLe &~D ermer 2007), the frac- 
tion of th e jet energy used to accele rate baryons (versus elec- 
trons; e.g. ISironi & Spitkovskvl 1201 ll) . and the slope of the in- 
jected UHECR spectrum. It is thus interesting to note that a shal- 
low injected energy spectrum both alleviates GRB energetic con- 
straints (j j4.2| l and may be necessary to fit the UHECR spectrum 
and composition measured by Auger if t he injected compositio n 
is indeed dom i nated by heavy nuclei (e.g. lHooper & Tavlorll2010h . 
iMetzger et al.l d201Ct) found that observed GRB spectra were best 
understood in the proto-magnetar model if magnetic reconnection 
( 33. U was responsible for powering the prompt gamma-ray emis- 
sion rather than shocks ( 33.2^ It is thus also important to note that 
numerical studies of magne tic reconnection indeed tend t o predict 
flat accelerated spectra (e.g. lRomanova & "Lo velace 1 992T) . 

Because the electron fraction in GRB outflows is uncertain, 
we cannot determine whether the heavy nuclei synthesized in GRB 
outflows are dominated by Fe group elements (A ~ 40 - 60) or 
whether the distribution extends to even heavier nuclei (A > 90). If 
the latter are present in at least a modest subset of events, a unique 
prediction of our model is that the UHECR composition may con- 
tinue to increase to nuclei heavier than Fe at yet higher energies 
> 10 20 eV. Making a measurement of the composition that is suf- 
ficiently accurate to test this prediction will, however, require both 
better statistics and a better understanding of the hadronic physics 
used to interpret the air showers. 

Another consequence of an ultra-heavy composition is that 
the accessible distance of sources may be appreciably larger than 
for Fe or protons (see 34.11 and eq. (8)). This alleviates energetic 
constraints on the GRB model for UHECRs, provided that the in- 
tergalactic magnetic field is not too strong ( 34.2l l. Although heavy 
nuclei probably compose a substantial fraction of the mass X h ~ 1 
in magnetized GRB outflows, the remainder is locked into 4 He. 
Helium is easily disintegrated into protons by the EBL and hence 
could contribute to the proton flux near the ankle. Alternatively, 
the proton-rich composition measured near the ankle could be sec- 
ondary particles produced by the disintegration of heavier nuclei 
by the EBL, or they could represent an entirely different source of 
UHECRs (Galactic or extra-galactic). 

If heavy nuclei from GRBs are indeed an important source 
of UHECRs, this would have several important consequences 
for GRB physics. For one, it would imply that GRBs outflows 
are magnetically-dominated, rather than thermally-driven fireballs 
(at least at their base). There is in fact growing evidence from 
Fermi observations that GRB jets may be magnetically-dominated 



(e.g. lZhang&Pe'erll2009h . If heavy nuclei are accelerated in GRB 
jets, this would also disfavor models in whic h GRBs are powere d 
by heating from neutron-proton collisions dBeloborodovl |2010|) . 
The high densities in magnetically-driven outflows make in un- 
likely that free nuclei will avoid being captured into heavy nuclei. 
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